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Abstract 

In type-II superconductors, the magnetic field enters in the form of vortices; their flow under application of a current 
introduces dissipation and thus destroys the defining property of a superconductor. Vortices get immobilized by 
pinning through material defects, thus resurrecting the supercurrent. In weak collective pinning, defects compete 
and only fluctuations in the defect density produce pinning. On the contrary, strong pins deform the lattice and 
induce metastabilities. Here, we focus on the crossover from weak- to strong bulk pinning, which is triggered either 
by increasing the strength f p of the defect potential or by decreasing the effective elasticity of the lattice (which 
is parametrized by the Labusch force /Lab). With an appropriate Landau expansion of the free energy we obtain 
a peak effect with a sharp rise in the critical current density j c ~ jo(ao£ 2n p)(£ 2 / a o)(/p//Lab — l) 2 . 
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Pinning of vortices by material defects is crucial 
in establishing the dissipation-free flow of a super- 
current. In recent years, the focus has been on weak 
collective pinning theory [1] describing the action 
of many competing defects; pinning then is due to 
fluctuations. The critical current is determined by 
the statistical summation of the competing pins [2] 
and is characterized by a quadratic dependence on 
the defect density. On the contrary, first attempts 
to describe vortex pinning go back to Labusch [3] , 
who studied the action of strong individual pins. 
Strong pins deform the lattice [4,5,6,7] and induce 
metastabilities; the pinning energy landscape be- 
comes multi- valued, producing a non-zero average 
of the pinning force and a critical current which is 
linear in the defect density. These two approaches 
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have been related in a recent work, where we have 
studied the pinning diagram exhibiting crossovers 
between various regimes when varying the defect 
density n p measured with respect to the vortex 
density 1/oq and the pinning force / p measured 
with respect to the elasticity, cf. Fig. 1. At low 
defect density, vortex lattice (bulk) pinning is rel- 
evant; for weak defect strength pinning is collec- 
tive, whereas defects with a force larger than the 
Labusch force f^ab (see below) pin the lattice indi- 
vidually. At higher defect densities we have estab- 
lished two regimes where pinning acts on individ- 
ual vortex lines. 

In this work, we focus on the crossover from 
weak- to strong bulk pinning, which is triggered ei- 
ther by increasing the defect strength k ~ / p /£ (£ 
is the coherence length of the superconductor) or 
by decreasing the effective elasticity C of the lat- 
tice; the critical defect force at the crossover is the 
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Labusch force /Lab = C£. The crossover is struc- 
turally related to the Landau theory of first order 
phase transitions and can be analyzed with an ap- 
propriate Landau expansion of the free energy. As 
a result, we obtain a strong increase in the critical 
current density at the crossover, 
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where jo is the depairing current density and ao = 
(Qo/B) 1 / 2 is the mean vortex spacing ($o is the 
superconducting flux unit); see below for the pre- 
cise definition of / p and /Lab- The crossover from 
weak- to strong pinning can be triggered by reduc- 
ing the effective elasticity C cx /Lab, producing a 
peak in the current density [1]. 

Below, we will describe the strong pinning situ- 
ation in more detail and use a Landau type of ex- 
pansion to find the jumps in the free energy land- 
scape. These jumps are then used to derive the re- 
sult (1) for the critical current density. 

In order to derive a quantitative criterion for 
strong pinning, we consider the effect of a single 
(strong) defect on the vortex lattice. We assume a 
defect located at the origin with pinning potential 
e p (r) producing the pinning contribution 

E p (r, u) = J2 e p (r)S 2 (R - R„ - u(R„, z j) 

to the free energy density of the vortex system; 
r = (R, z) and the vortices are positioned at R„ — 
u(R v ,z), with R„ the equilibrium positions and 
u the displacement field. The clastic part of the 
vortex lattice free energy reads [2] 
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l (-k) [G a0 (k)] V(k), (2) 



with G af3 (k) the Fourier transform of the elastic 
Green function G"^(r); Greek indices denote the 
in-plane components (x, y) and we sum over double 
indices. 

The displacement field u is obtained from vari- 
ation of the total energy (elastic and pinning) , 
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Fig. 1. Pinning diagram delineating the crossovers between 
various pinning regimes involving collective versus individ- 
ual pinning and ID-line versus 3D-bulk pinning: 3D wep 
- bulk weak collective pinning, ID cp — collective line pin- 
ning, ID sp - strong line pinning, 3D sp - bulk strong 
pinning. The crossover from weak- to strong bulk pinning 
occurs at the Labusch force /Lab- 

with r„ = (R„, z) and f p = — V u e p (u) the pinning 
force of the defect. In the last equation, we have 
assumed a defect of range much smaller than a 
pinning at most one vortex; we further have cho- 
sen Rd as the distance to the vortex closest to the 
defect. Evaluating (3) for r„ = (Rd,0), we arrive 
at the self-consistency equation 

u a (R,0) = G af} (r = 0)fg\R + u(R,0),0] 

-C- 1 / p Q [R + u(R,0),0] (4) 

with the effective elastic constant 
d 3 k 
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For ao < A the effective elasticity is C ~ so/ao 
with £o = ( < I ) o/47rA) 2 the vortex line energy. The 
solution at (R, 0) allows for the calculation of the 
whole displacement field u(R [/ , z), 



u(R», z) = G afj (R v - R, z)Cu(R, 0), 



(6) 



and we can calculate the total free energy 
e p in(u, R) (containing both, the contribution from 
elasticity and disorder) of the vortex system as a 
function of u(R, 0), 
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= G a!i (R„ - R d , z)f% [R d + u(R d , 0) , 0] , 
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Fig. 2. The pinning energy e p ; n and the force / p j n are 
sketched as a function of the drag Rd between the vortex 
system and the defect. In the case of weak pinning (dashed 
line), the two quantities are single- valued (hence the aver- 
age over the pinning force is zero), whereas in the strong 
pinning case (solid line, the dotted parts denote unstable 
branches) the average over the force is determined by the 
jump Ae p ; n in the pinning energy. The thick lines show 
the choice of the two possible branches that occurs when 
dragging the lattice through the defect along the positive 
x direction. 



^Cu 2 + e p (R + u). 



(7) 



For weak pinning, the displacement u is small 
and the solution u « f p (R)/(7 of (4) is unique; 
the pinning energy e p \ n (R) = e p ; n (u(R),R) then 
is a single-valued function. Strong pinning, how- 
ever, produces multi-valued solutions u(R, 0) and 
e p in(-R), as can be seen in Fig. 2. The displacement 
turns multi- valued when 8ru — > oo. Assuming a 
defect symmetric in the plane (e p (R) = e p (R)) 
and dragging the lattice through the defect center 
along the x-axis (Rd = (x, 0)), we find 



d x u{x) 



C-f'(x + u) 



(note that x > implies u < 0). The displacement 
collapses when the (negative) curvature / p of the 
defect potential overcompensates the elasticity of 
the lattice, 



This (Labusch) criterion [3] for strong pinning in- 
volves the maximal negative curvature above the 
inflection point; it tests an individual isolated pin 
and classifies it as weak or strong. 

In a next step, we calculate the average pin- 
ning force in the strong pinning case. Using the 
self-consistency equation (4), we calculate the 
derivative of the pinning energy e p i n (x,y) = 
e P in(u(R), R) with respect to the drag parameter 
x, 

d x e pin (x, y) = Cuc^u - f p (R + u)[x + d x u] 

= -/*(R + u); (9) 

this equation relates the force along x exerted by 
the pin at the vortex position R + u(R, 0) to the 
free energy landscape containing both, the contri- 
bution from pinning and elasticity. The pinning 
force has to be averaged over defect locations and, 
due to the bistabilities, depends on the prepara- 
tion of the system; here, we focus on the critical 
current density and therefore search for the maxi- 
mal force against drag. Averaging the accumulated 
drag force over the 'impact parameter' y, we obtain 
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the integral over x can equally well be interpreted 
as an integrated drag force or as an average over pin 
locations (maximized with respect to the bistable 
branches). We express the average along the x-axis 
through the jump Ae p i n (y) in the pinning energy, 
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(ii) 



d x fp(x + u) = -d 2 x e p {x + u) = C. 



(8) 



where we have assumed a maximal trapping dis- 
tance t± along the y axis and require the pin not to 
overdrag the vortex; then the periodicity of e p i n is 
the same as the lattice periodicity and a = a . The 
remaining task is the determination of the jump 
Ae p i n ; as we will see below, it does not depend on 
the impact parameter y. 
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Fig. 3. Displacement u and energy e p j n versus drag x of the 
vortex lattice relative to a defect producing a Lorcntzian 
pinning potential e p L = — e ()/(l + " 2 /5 2 ) ( wc measure 
units of length in £ and units of energy in eo ) . The evolution 
of the bistability appearing above the Labusch criterium 
is shown for pinning centers of increasing strength with 
C'/k = 0.85,0.75,0.65,0.55; dashed lines denote unstable 
branches. The inset shows the geometry defining the drag 
parameter x and the displacement u. 

In a first step, we discuss the solutions of the self- 
consistency equation (4) when dragging the lattice 
through the defect center along the positive x di- 
rection and calculate the displacement u numeri- 
cally for a Lorentzian defect potential 

The result is shown in Fig. 3, where we have plot- 
ted the displacement u and the pinning energy e p ; n 
for different values of C/k ~ /Lab// P , where —k = 
— eo/2£ 2 is the maximal negative curvature of the 
defect potential (realized at u = u K = £), which 
is the relevant curvature in the Labusch criterion 
(8). The branch u_ (u + ) corresponds to the pinned 
(unpinned) vortex and becomes unstable at x_ 
(x + ); the dashed lines denote unstable branches. 

The solutions u(x) of the self-consistency equa- 
tion (4) correspond to the minima of the functional 
e p i n [w, a;] for fixed x, cf. Fig. 4. Starting far from 
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Fig. 4. Schematic plot of the pinning energy as a function 
of the displacement u for different (negative) values of 
the drag x. The sketch shows the trapping process of an 
unpinned vortex at the point of instability —x+ and the 
appearance of the jump Ae p ; n . 

the defect with a negative drag parameter x, there 
is only one solution u + corresponding to the un- 
pinned vortex system, cf. Figs. 2 and 4. With the 
vortex approaching the defect a second minimum 
appears, but is separated from the first one by a 
barrier. The relevant solution therefore remains the 
unpinned one, although at some point the pinned 
solution becomes lower in energy. For x = —x + the 
unpinned solution becomes unstable and the vor- 
tex gets trapped, as shown in Fig. 4. It remains so 
until at x = x_ the pinned solution u_ turns un- 
stable and the defect releases the vortex. The cor- 
responding jumps in the free energy are shown in 
Fig. 2. 

In order to attain a more quantitative descrip- 
tion of the different branches and their correspond- 
ing jumps, we perform a Landau type expansion 
of the free energy close to the Labusch condition 
k = C and use the analogy to the Landau theory of 
first order phase transitions. We denote the loca- 
tion of the maximal negative curvature —nhyu K 
and expand the pinning potential e p (x + u) around 
this point, 

e p (x + u) w — e + v(x + u — u K ) (12) 
--(x + u-u K ) +—{x + u-u K ). 

For the Lorentzian potential the expansion param- 
eters read u K = £, e = e /2, v = e /2£, k = e /2£ 2 , 
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and a = 3eo/£ 4 . We insert this expansion into the 
pinning energy e p j n = Cu 2 /2 + e p (R + u); the lat- 
ter maps to the free energy of a one-component 
(Ising) magnet in a magnetic field [8] if we define 
the order parameter <j> = x + u — u K , the reduced 
'temperature' r = C — n, and the 'magnetic field' 
h = C(x — u K — v/C), 

e pin [u,x] = e mag [</>, h] = ^ 2 + ^-</» 4 - h(j), (13) 

where we have dropped the constant term. Con- 
sider the zero field case first: In the paramagnetic 
phase at high-temperatures r > the order pa- 
rameter <j) vanishes, while <f> oc d= |-r | x / 2 in the 
low temperature (r < 0) ferromagnetic phase. 
The two ferromagnetic states (</> > and <j> < 0) 
are separated by an energy barrier increasing as 
Ae mag ocr 2 . In the low temperature phase r < 0, 
a field h induces a first order transition across 
h = between the two ferromagnetic states. 
Metastable/hysteretic behavior appears within 
the field regime \h\ < h* ', with h* oc |t| 3 / 2 follow- 
ing from comparing the barrier Ae mag and the 
field energy h(f>. 

The above considerations translate to the pin- 
ning problem in the following way: the high 'tem- 
perature' phase t > describes weak pinning. The 
two low temperature states stand for the pinned 
(<f> < 0) and unpinned (4> > 0) configurations, 
which transform into one another via the first or- 
der transition. The hysteretic regime translates to 
the bistable pinning domain bounded by 

x± = u K + {v T- h*)/C, 

with h* = (2/3C") v /27"|t| 3/2 - At the instability 
x = x- (h = h*) the order parameter jumps 
from the pinned solution u_ = (f>- + u K — X- to 
the unpinned state u+ = 4>+ + u K — x_, where 
(f>- - -y^H 1 / 2 and 4>+ = 2 v /27^|r| 1 / 2 . As 
a result, we find the jump Ae = e p i n (u_,a;_) — 
e P i n (w+, X-) = (9/2a) r 2 in the free energy; a jump 
of equal magnitude is found at x+. The total jump 
Ae p i„ in the pinning energy then reads 

Ae pin = 2Ae = (9/a) (C - k) 2 . 

Within this Landau description, the Labusch con- 
dition C = k corresponds to the critical end-point 
T = T C terminating a line of first order transitions. 



Going beyond this one-dimensional analysis, we 
have to account for a finite 'impact parameter' y of 
the vortex with respect to the defect (Rd = (x, y)) 
and determine the transverse trapping distance t± 
for the vortex to get trapped when dragged along 
x, cf. Fig. 5. Making use of the rotational symme- 
try of the problem, we find circles with radius x± 
limiting the bistable regions. The jumps then occur 
when a solution is dragged through its instability 
line: The unpinned vortex gets pinned when x = 
— (x 2 ^— j/ 2 ) 1 / 2 ; this only happens for impact param- 
eters \y\ < x+, which leads to a trapping distance 
tj_ = 2x + (note, that t± w u K + v/C remains finite 
at the transition to weak pinning). The pinned vor- 
tex is released when x = {x 2 _ — y 2 ) 1 / 2 . Note, that 
the corresponding magnitudes of the jumps do not 
depend on y because of the planar symmetry. 

| uncompensated 
y pinned 




vortex trajectory 



Fig. 5. Trapping area around the defect for a circular sym- 
metric situation with pinned, unpinned, and bistable re- 
gions; the dark uncompensated area produces the net pin- 
ning force. A vortex trajectory for finite impact parameter 
y < is shown. 

We insert the value for the jump Ae p ; n into (11), 
use t± = 2x + and obtain the pinning force (/ p i n ) = 
— 2x + Ae p i n /ag. Adding the contributions of inde- 
pendent strong defects with density n p , the total 
pinning force density reads F p i n = n p (/ p ; n ) and 
the critical current density is obtained as j c ~ 
—cF p i n /B. Close to the Labusch condition, we fi- 
nally obtain the critical current density 

Jc- Jo a o en p ^ (t^-i) , (14) 

% V/Lab / 
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where we have used a defect size £, implying t± ~ 
£; we define the relevant defect strength / p = £k = 
£ max„ [/p(u)] in terms of the maximal curvature of 
the defect potential. The Labusch force then takes 
the form /i^b = £C ~ £o£/ a o- Interpolating with 
the weak collective pinning result [2] 




(we assume pinned bundles of size larger than 
A), we observe a sharp rise in the critical current 
density once the strong pinning force overcomes 
the weak pinning result; this peak effect [1] oc- 
curs close to the Labusch condition at / P //Lab = 
1 + {ao/\)y/ao£ 2 n p . 

The crossover from weak- to strong pinning can 
be realized in experiments: increasing the magnetic 
field and approaching the upper critical field H C2 
leads to a marked softening of the elastic moduli, 
thereby reducing /i^b oc C; the resulting crossover 
into the strong pinning regime then is observed as 
a peak [1,9] in the current density. Whether such a 
crossover from weak- to strong pinning can explain 
the recently observed peak effects (see, e.g., [10,11] 
and Refs. therein) remains to be seen. 
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tional Foundation. 
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